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Abstract 

It is shown that if the decoherence matrix corresponding to a qubit 
master equation has a block-diagonal real part, then the evolution is de- 
termined by a one-dimensional oscillator equation. Further, when the full 
decoherence matrix is block-diagonal, then the necessary and sufficient 
conditions for completely positive evolution may be formulated in terms 
of the oscillator Hamiltonian or Lagrangian. When the solution of the os- 
cillator equation is not known, an explicit sufficient condition for complete 
positivity can still be obtained, based on a Hamiltonian/Lagrangian in- 
equality. A rotational form-invariance property is used to characterise the 
evolution via a single first-order nonlinear differential equation, enabling 
some further exact results to be obtained. A class of master equations is 
identified for which complete positivity reduces to the simpler condition 
of positivity. 

1 Introduction 

Master equations are useful for representing the evolution of non-isolated quan- 
tum systems, where the details of the interaction with the environment are 
encoded in a decoherence matrix or memory kernel pQ. In practice, given the 
complicated nature of typical environments, one must often work with master 
equations that have been derived via approximations and/or phenomcnologi- 
cal considerations. Unfortunately, in such cases the resulting master equations 
may not correspond to evolutions that are physically possible - for example, the 
density operator may evolve to have negative eigenvalues. 

The distinction between physical and nonphysical master equations is not 
obvious in general. For example, Diimcke and Spohn have pointed out that 
subtly different ways of approximating the physical principle of weak coupling, 
for a qubit interacting with a thermal reservoir, can variously lead to either 
physical or nonphysical evolution |2J. Similarly, Barnett and Stenholm have 
shown that the assumption of an apparently innocuous exponential memory 
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kernel, describing a harmonic oscillator coupled to a reservoir, leads to negative 
probabilities [3]. 

There is, therefore, interest in finding conditions on master equations which 
ensure that the corresponding evolution of the density operator is physical 
[TJ IH O [HI M QZ3- In this regard, it is not sufficient to merely ensure 
that the density operator remains positive under evolution. In particular, if 
some auxiliary system is correlated with the system of interest, but does not 
interact with it, then the corresponding joint density operator must also remain 
positive under the evolution. This requirement is stronger than positivity, due 
to the remarkable nature of quantum correlations, and is called complete posi- 
tivity [TJ llll 112] . Determining conditions for qubit master equations to generate 
completely positive evolution is the focus of this paper. 

The general form of a memoryless master equation, for a qubit system de- 
scribed by density operator p, is (eg, equation (2.7) of [4]) 

p = A t (p) := -i[H(t),p] + (l/4)2^7ife(t) {2a 2 pa k - a k ajp - pa k aj) . (1) 

Here H (t) is a Hamiltonian operator, the Uj are the Pauli spin matrices, and 
7(i) is a 3 x 3 Hcrmitian matrix which will be referred to as the decoherence 
matrix. Master equations written in memory-kernel form can also be reduced 
to the above time-local form, provided that a particular inverse exists |T3] • 
Note that the (minimal) Lindblad form of the master equation corresponds to 
the eigenvalue decomposition 7 = J2i=i A/e^e^t of 7, i.e., one has dHH[S] 
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p = -i[H, p] + (1/4) A/ (2 Uph\ - L\Up - ph\h^ , 
\=\ 

with Li := • a (thus, tr[i]L fc ] = 2 8 jk ). 

In the case of no explicit time-dependence, the necessary and sufficient con- 
dition for completely positive evolution is simply that the decoherence matrix 
is nonnegative [4, 5 , i.e., 7 > 0. While this remains a sufficient condition when 
H and/or 7 depend on time (since the evolution is then a composition of a 
sequence of infinitesimal completely positive evolutions), finding the necessary 
and sufficient conditions for complete positivity in the time-dependent case is a 
very difficult problem. 

This problem has been solved by Wonderen and Lendi 6 , and independently 
by Maniscalco [8), for the case 

/ 7 ig \ 
H(t) = (l/2)ha 3 , 7 (i)= I -ig 7 , 

V 73/ 

which is applicable to several systems of physical interest (including an example 
where the master equation is obtained from a memory- kernel form [8]). The 
more trivial case of characterising complete positivity when the decoherence 
matrix is diagonal and H — is also solvable, and is reviewed in [7j. 
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In this paper the more general form 

/ 7i f + ig ir \ 

H(t) = (l/2)ha 3 , 7 (*)= f-ig 72 is (2) 

\ -ir -is 73 / 

is considered, where /, g, h, r, s and the 7j are all real functions of time. Thus 
7n 7^ 722 in general; 712 may have a real component; and (more trivially) there 
is no restriction on the imaginary part of the decoherence matrix. Note that 
this form is equivalent to the condition that the real part of the decoherence 
matrix is block-diagonal. As will be seen in section 2, it is also equivalent 
to the condition that the damping matrix in the Bloch representation is block- 
diagonal, i.e., equation ([2]) corresponds to the case that damping in one direction 
is decoupled from damping in the remaining two orthogonal directions. 

It is shown in section 3 that the evolution corresponding to equation ^ 
is in general determined by the solutions of a one-dimensional time-dependent 
oscillator equation 

d 2 q/dT 2 + k{r)q, (3) 

where the reparameterised time r and the 'spring constant' k(r) are determined 
by H and 7. Explicit solutions for the oscillator motion yield explicit solutions 
for the corresponding master equation. 

The general form in equation ([2]) satisfies the property that the master equa- 
tion is form-invariant under arbitary (time-dependent) rotations of the system 
about the z-axis. This provides a useful gauge-like degree of freedom for sim- 
plifying the master equation, which is exploited in section 4 to characterise 
the qubit evolution via a single nonlinear first-order differential equation. This 
equation may be solved, for example, when 

h = 0, f = Kty - 72) 

in equation ([2]), for some constant K, which generalises the abovementioncd 
previously considered cases El [8] , and includes the case of a two-level atom 
in a time-dependent squeezed vacuum. 

In section 5 it is shown that when the full decoherence matrix is block- 
diagonal, i.e., when 

r = s = (4) 

in equation |(5J) (corresponding to any 'drift' in the Bloch representation being 
confined to the z-direction), then the necessary and sufficient conditions for 
completely positive evolution may be formulated in terms of either the oscillator 
Hamiltonian or Lagrangian, depending on whether k is positive or negative. 
This leads to an explicit characterisation of complete postivity, whenever the 
oscillator system is solvable. 

Further, Lagrangian and Hamiltonian oscillator inequalities are derived which 
lead to a useful sufficient condition for complete positivity. In particular, this 
condition is formulated solely in terms of the master equation parameters, and 
hence can be applied whether or not the corresponding solution is known. It is 
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quite strong, being both necessary and sufficient in a number of cases, and is 
invariant under rotations of the system about the z-axis. 

Finally, a class of master equations is identified for which complete positivity 
is equivalent to positivity, corresponding to the case of unital evolution with 
no damping in the z-direction. This class includes a group of Bloch-Redfield 
master equations in a white noise limit, recently considered by Whitney which 
generate positive evolution [9]. It immediately follows that the evolution is in 
fact completely positive for this group. This case is of further interest in that 
the question of complete positivity can be settled despite being unable to solve 
for the evolution explicitly. 



2 Evolution in the Bloch representation 

It is convenient to rewrite the master equation |T]) in terms of the Bloch vector 
v, where 

p=(l/2)[l+v-4 (5) 

Substituting into the first equality in |T]) and taking the trace with Oj then leads 
to the equivalent Bloch equation 

v = u + Dv, 

where 

Uj := (l/2)trfoA(l)], D jk := (1/2) tr[<7,-A(<7 fe )] 

arc termed the drift vector and damping matrix respectively [8] . This is a first- 
order inhomogenous differential equation, and hence the general solution is of 
the form 

v(t) =Af(t)v(0) + w(t), (6) 

for some matrix M and vector w. Substitution of © into the Bloch equation 
yields the equivalent evolution equations 

M = DM, w = u + L>w, 

for M and w, subject to the inital conditions 

M(0) = I, w(0) = 0. 

Now, to first-order in e, 
det M (t + e) = det M det (J + eD) = dct M + e£>#) = det M(l + etr[D]), 

3 

and hence it follows that 



det M = exp (J dstr[D(s)] \ > 0. 
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Thus the inverse matrix M 1 always exists, and it is straightforward to check 
that the solution for w is given by 



w{t)=M(t) dsM- 1 (s)u(s). 



(7) 



Hence, solving the general master equation ([I} is equivalent to solving the matrix 
equation 

M = DM, M(0) = / (8) 

for the evolution matrix M . Similar considerations apply to higher-dimensional 
quantum systems. 

For the general qubit master equation in |T]), the above definitions of u and 
D and the relations 

tr[a a a b a c ] = 2ie abc , tv[a a a b a c a d ] = 2 {5 ab 5 cd + 5 ad 5 bc - 5 ac 5 bd ) (9) 

may be used to calculate 

ui = -Im{7 23 } et cyclic, D jk = Re{7 jfc } - S jk tr[j] - ^ e jkl tr[Hai\. 

i 

It follows that, for the particular forms of H and 7 defined in ^j, the drift 
vector and damping matrix are given by 









"I 




-I 









-72 - 73 f-h 
f + h -71 - 73 
-71 - 72 



(10) 



It is the block-diagonal form of the damping matrix D, corresponding to damp- 
ing in the z-direction being decoupled from damping in the x and y directions, 
that provides the basis for the main results of this paper. 



3 Reduction to an oscillator system 
3.1 Oscillator form 

From equations JHJ) and (TT0|) . the evolution matrix M is itself block-diagonal, 
i.e., one has 

(xi x 2 \ 
yi V2 . (11) 
A J 

Substitution into the evolution equation ([5]) yields in particular that 

dA/dt = -( 7l + 72 )A, dA/dt = -( 7l + 72 + 2 73 ) A, 
where A denotes the subdeterminant x\yi — x^yi, and thus 



A = exp 



ds (71 + 72) 



A = A exp 



-2 / ds 73 
Jo 



(12) 
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Hence only three independent parameters of M remain to be determined. 
Consider now the 2-vector equation 



-72 - 73 

f+h 



f-h 

-71 - 73 



where (x, y) — (xj,yj) (j — 1, 2), corresponding to the evolution of the upper 
block of M. Defining the quantities 

q:=xeIo ds ^ 2+la \ p := X ds ^^y, G := X ^ , (13) 

one immediately finds that 

q=(f-h)G~ 1 p, p=(f + h)Gq, 
and hence the evolution of q and p is described by the quadratic Hamiltonian 

H(q,p,t) -=\{f- h)G- V -\{f + h)Gq 2 . 

To obtain the canonical oscillator form, note that the corresponding action, 
J dt (pq — H) , may be rewritten as 



dt [(/ -hy'Gq 2 + {f + h)Gq 2 



dr 



kq 2 



providing that one defines the functions r(t) and fc(r) via 

t:=(f-h)G-\ r(Q)=0, k := G 2 (h + f)/(h - /). 



(14) 



Hence, with respect to the reparameterised time r, one has the oscillator equa- 
tion 



d 2 q/dr 2 + kq = } 



(15) 



as previewed in the introduction. 

It follows that when the oscillator equation (fT5| can be solved, then the 
evolution matrix M can be determined. In particular, such a solution must link 
the oscillator state to its inital state via a linear relation of the form 



q 

dq/dr 



a b 

c d 



(dq/dr)o 



(16) 



Since dq/dr = q/i = p from the defining equations (TT3")) and (Ti~4"|) . it follows 
further via (H]) and HID that 



- J ds (72+73) Q 

e~ / ds ( 7l+73 ) 



a b 
c d 



Vq 



G 



Recalling the initial condition M(0) = /, the corresponding explicit form 



M 



I ae- $ ds{l2+Ti) 
ce ~J ds (71+73) 

V 



be~f ds (72+73) 




(17) 



■J ds (71+72) 



is obtained for the evolution matrix. Thus, the master equation is solvable 
whenever the oscillator matrix in (1161) is known. 



3.2 Examples 

As a simple example, consider the case / = = 71 — 72. Then, k = G = 1 
and r = — J Q ds h(s) from equations (fl3|) and (| 14[> . and hence from (fT5|) the 
oscillator matrix is 



6 \ _ / cost sinr \ _ I cos f*dsh(s) —sin f*dsh(s) 
d J ~ \ -sinr cost / I s i n / 'ds/i(s) cos f*dsh(s) 



The solution of the corresponding master equation in the Bloch representation 
then follows via equations 10 -0 and ([17]). Note this case corresponds to that 
considered previously by Wonderen and Lendi [5] and Maniscalco [8], It is 
significantly generalised in the following section. 

As a second example, consider the case h = = 71 — 72, corresponding to a 
symmetric damping matrix in the interaction picture. One finds k = — 1, G = 1 
and t = J Q ds f(s) and hence that 




cosh J^ds f(s) sinh J* Q ds f(i 
sinh J Q ds f(s) cosh J Q ds f(i 




Finally, it proves useful to consider the degenerate case / = —h, for which 
one has k = and r = 2 JL ds/G _1 . Thus the corresponding oscillator system 
degenerates to free particle motion, with zero frequency, and 

U d)-{o i)-[o h 1 J ' 

Note that this reduces to the identity matrix when one further has / = 0, 
corresponding to the trivial case of a diagonal damping matrix in equation pop 
(previously reviewed in [7]). More significantly, this solution is used in the 
following section, together with rotational form-invariance, to further reduce 
the evolution of the master equation to a single nonlinear first-order equation. 
Note for all the above examples that 

ad-bc=l. (18) 

This holds more generally, as may be derived directly from the oscillator equa- 
tion (TH)]) . or via equations (fT2")) and (TTTj) . Note also that the oscillator equation 
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breaks down for the singular case / = h, since k becomes undefined in equation 
(fT4|) . However, one can solve equation ([8]) for M directly in this case, to find 
it corresponds to an oscillator matrix equal to the transpose of the one given 
above for / = —h. 



4 Exploiting rotational form- invar iance 
4.1 Rotational form-invariance 

Consider now a rotation of the system about the z-axis, by some time-dependent 
angle a{t). Such a rotation corresponds to the unitary transformation p :— 
U pW of the density operator, with U := exp[—ia(t)cr 3 }. Denoting the corre- 
sponding rotation matrix by R, one has U'crjU — J^ fc Rjk&k, and substitution 
into (JlJ yields the transformed master equation 

p = At(p) := -i[H(t),p] + (l/4)^7,-k(t) (2<Jjpa k - a k afp- ~pa k u 2 ) , 



with 



H = H + aa 3 , x f = R 1 R T . 




A tilde will be used in general to denote quantities in the rotated frame. 

For the particular class of 'Bloch-diagonal' master equations described by 
equation (|2|), one finds in particular that 



H = (l/2)ha 3 , 7=| / -.. | . 

where 

71 := 71 cos 2 a + 72 sin 2 a — f sin 2a, 72 := 71 sin 2 a + 72 cos 2 a + f sin 2a, 

h = h + 2d, f = / cos 2a + (1/2) (71 -72) sin 2a, 

f — r cos a — ssina, s — r sin a + s cos a. 

Note that g, 71 -I- 72 and 73 do not change under the rotation. 

Comparing equations ^) and (fTO|) it is seen that the form of equation © 
is preserved by such rotations. Moreover, from equation and the property 
Wall — Ra, the Bloch vector transforms under such rotations as 

v = Rv. 

Substitution into equation ^ then yields the relations 

M = R T AIR(0), w = R T w. (20) 

Thus, if the master equation can be solved in the rotated system, by a judicious 
simplifying choice of the function a(t), then the solution with respect to the 
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original system can also be determined. This rotational degree of freedom allows 
non-trival new exact solutions to be obtained for qubit master equations, via 
reduction of the evolution to a single first-order differential equation for a, as 
is shown below. It is also relevant to the discussion of complete positivity in 
section 5, where it is used to obtain a rotationally-invariant sufficient condition. 



4.2 Reduction to a first-order nonlinear equation 

For any master equation of the form define the 'zero frequency' gauge or 
picture via the condition 

/ = -h. 

Note that this is rather different from the standard 'interaction' picture, which 
corresponds to h = 0. From the above expressions for / and h, this condition 
may be rewritten as the first order differential equation 

a+-[2f cos 2a + ( 7l - 72 ) sin 2a + 2h] = (21) 

for a. 

Now, the choice f = —h corresponds to the degenerate case of a zero- 
frequency oscillator considered in section 3.2, and therefore the corresponding 
oscillator matrix can be immediately written down as 

(a b\ = ( 1 2 /* ds f Jo du } 

\ c d J y " 1 

Hence, if a can be determined from equation (|2ip . then M can be determined 
via (|17[) . and the evolution matrix M for the original master equation follows 
via (|20|) . 

Thus, remarkably, solving the master equation is equivalent to solving the 
zero-frequency gauge equation (|2ip for a. Naturally enough, solving this equa- 
tion explicitly cannot be done in general, as it would amount to solving a general 
time-dependent oscillator problem. 




4.3 Example: a new solution 

It is possible to solve equation (|2"Tj) in some cases of interest other than the 
examples of the previous section. For example, consider the case 

h = 0, / = tf( 7l - 72 ) (23) 

for some constant K. This case corresponds to fc(r) being explicitly time- 
dependent in the oscillator equation (|15p . and includes the known cases 7 i — 7 2 = 
/ = h = [6l[S and / = h = [7] in particular. However, it is rather more gen- 
eral, including, for example, a two-level atom with natural linewidth 7 coupled 
to a squeezed vacuum described by squeezing parameter = r(t)e l9ot , with 
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(slowly varying) time-dependent squeezing parameter r(t) and fixed squeezing 
angle 6q. This corresponds to the particular choice [1] 

7i + 72 = 7 cosh 2r, 71 — 72 = — 7 cos #0 sinh 2r, / = —(1/2)7 sm ^0 sinh2r 

in equation with all other coefficients vanishing, and hence to K = — ^cot 9q 
in equation (]23[) above. 

To solve the master equation in this case, define the (constant) angle by 



cos20:=2ii:/Vl + 4is: 2 , sin 20:= 1/^/1 + 4K 2 . 
The 'zero frequency' equation ([21)1 can then be rewritten in the separable form 



(sin 20) j t 2(a - 0) = -(l/2)( 7l - 72) cos2(a - 0), 

which may be immediately be integrated to give, assuming that a(0) = for 
convenience, 

(sin20)logtan(a - + tt/4) = -(1/2) / ds (71-72). 







Inverting gives the explicit expression 



a = — 7r/4 + tan 1 



o -(l/4,K)VT+4Kl /* (71-72) 



(24) 



for the gauge function a, as desired. Thus, for all master equations satisfying 
condition (|23[) , the explicit evolution in the Bloch representation can be obtained 
by (i) applying equation (|17p to equation (|22p for the above choice of a, to obtain 
M , and (ii) finding the evolution matrix M via relation ([2"D"[) . 

5 Complete positivity for block-diagonal 7 

5.1 Necessary and sufficient conditions 

The evolution of a quantum system, described by some linear map pt — <j>t(po), 
is completely positive if and only if the corresponding Choi matrix C has no 
negative eigenvalues, i.e., if and only if C > [14]. For the case of qubits, 
it is convenient to calculate this matrix C with respect to the basis set used 
in section 4 of reference [7] (corresponding to the matrix S^ w ^ therein), and 
multiply by a factor of 2, so that the Choi matrix is the 4x4 matrix defined by 

C jk := (1/2) F m 

m, n 

where F mn := (l/2)tr[a m 0(tr n )], and the indices run over 0, 1, 2, 3 with ao := 1. 
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From the Bloch representation in equations ([5]) and ©, the evolution map 
4> is given by 

(j>{X) = (l/2)tr[X] (1 + w • a) + (1/2) (Mtv[Xa\) ■ a 

where a denotes the 3-vector (<7i, <j 2 , 03), and hence Foo = 1, Fqj — SjQ, Fjo = 
wj, and Fjk = Mjk for j,k = 1,2, 3. It follows, using properties © of the Pauli 
matrices, that the coefficients of the Choi matrix are given by 

Coo = 1 + tr [M] , C 01 =w 1 +i (M 23 - M 32 ) , 

Cn = 1 + Mn - M 2 2 - M33, C12 = M12 + M 2 i + iw 3 , 

with the remaining coefficients determined via cylic permutations of 1,2,3 and 
C = . Checking positivity of the Choi matrix for a general qubit evolution 
essentially requires finding the singular values of M [15] . and hence there is no 
general explicit condition in terms of the coefficients of M and w. 

However, for master equations with a block-diagonal decoherence matrix, the 
Choi matrix has a relatively simple form. This case is equivalent to r = s = 
in equation (|2|) and implies that, in addition to the damping matrix D being 
diagonal, the 'drift' vector u in (flQ|) is confined to the z-direction. Equations 
0, dTUJ) and (HID then lead to 

/ l + xi+y 2 + A w 3 + i{x 2 ~yi) \ 

q_ 1 + x-l - y 2 - A X2 + Vi + iw 3 

x 2 + yi - iw 3 l — xi + yz — A 

\ w 3 -i(x 2 -yi) l-Xi-y 2 + A J 

(25) 

for the corresponding Choi matrix, with 

w 3 = -A [ dsgA- 1 . (26) 

The condition C > thus reduces to the positivity of the two 2x2 sub- 
matrices composing C, i.e., to the positivity of the traces and determinants of 
these submatrices. This yields, after some rearrangement, the necessary and 
sufficient conditions 

A<1, S := x\ + xj + y\ + y% < 1 + A 2 - w\ - 2\A- A| (27) 

for complete positivity. Note from (TT2")) that A and A are explicitly defined in 
terms of the master equation parameters, as is w 3 (given above). 

The question of complete positivity therefore reduces to knowledge about 
the quantity S on the left hand side of second inequality in (|27[) . This quantity 
may of course be calculated when the solution of the master equation is known, 
such as for the examples in sections 3 and 4, thus completely determining the 
conditions for complete positivity in these cases. More generally, however, only 
partial conditions can be explicitly determined in terms of the master equation 
parameters, as discussed below. 
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5.2 Necessary conditions 



Here two necessary conditions for complete positivity are noted, for master 
equations having a block-diagonal dccohcrcncc matrix, which do not require the 
solution of the master equation. Both conditions are formulated in terms of 
quantities that are invariant under rotations about the z-axis. 
First, the condition A < 1 in (|27|) reduces via equation (| 1 2[) to 



ds (71 + 72) > 0. 



Second, noting that the quantity S in (|27[) is the sum of the squares of the 
singular values of the upper block of M in (fTTj) . and that the positive quantity A 
is their product, it follows via si + s 2 > 2siS2 and (|12[) that complete positivity 
requires 



A 2 -2A 



\l-e Jo "| + e Jo ' 3 



1 > 0. 



Note that this quadratic condition is certainly satisfied when the corresponding 
discriminant is negative, i.e., when 



I-. -2 f ds73i . —2 f* ds 73 . ~~ 

|1 — e J o I + e ' < J1 + 
which in turn is guaranteed when 

ft 

dsjz > 0. 



5.3 Sufficient condition from a Lagrangian inequality 

To obtain a nontrivial sufficient condition for complete positivity, it is convenient 
to begin by working in the interaction picture, so that h — 0. In this case one 
has k = —G 2 for the corresponding (inverted) oscillator system in section 3, and 
hence the oscillator Lagrangian is given by 

L{q,dq/dT,T) = (1/2) [(dq/dr) 2 + G 2 q 2 ] . 

Noting that the conjugate momentum is p = dq/dr, the value of the Lagrangian 
at any given time follows from equation (|16p as 

L = (l/2)G 2 (aq Q + b Po ) 2 + (l/2)(azo + d Po ) 2 . 

Hence, if L\ and L2 refer to the values of L at time t for the canonical initial 
states (qo,po) — (1,0) and {qo,po) = (0,1) respectively (actually, any two or- 
thogonal initial states of equal norm will do), then their average value evolves 
as 

I := (l/2)(Li + La) = (1/2) [G 2 (a 2 + b 2 ) + (c 2 + d 2 )] . 
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Note that L is invariant under phase space rotations. 

Comparing this expression with equations and (|17|) . and noting the 
definition of G in equation (fT3"|) , it follows that 

a 2,2,2,2 r»~F ~ 2 f* rfs (71+73) /r,a\ 

S = x l +x 2 + y 1 + y 2 = 2Le J° >. (28) 



Hence the complete positivity condition \2T\j may be interpreted as an upper 
bound on the average Lagrangian value of the corresponding oscillator system. 

In particular, any upper bound for L immediately generates a sufficient 
condition for complete positivity. One such bound is obtained here, using a 
generalisation of the method given by Boonserm and Visser for obtaining bounds 
for a 2 + b 2 + c 2 + d 2 (rather than for L) , relevant to one-dimensional scattering 
coefficients 16 . In section 5.3 this bound is shown to in fact be applicable to the 
case of arbitrary h, as a consequence of rotational invariance. This further allows 
a Hamiltonian upper bound to be obtained for the time-dependent harmonic 
oscillator. 

First, define the quantities X±, Z by 

X± := G(a 2 + b 2 ) ± G _1 (c 2 + rf 2 ), Z :=ac + bd. 
Note from the determinant property (fT5|) that 

X 2 - X 2 = 4(a 2 + b 2 ){c 2 + d 2 ) = A(ac + bd) 2 + 4(ad - be) 2 = AZ 2 + 4. 
Now, the oscillator equations (|15p and (fi"6| imply that 



d_ f a b \ _ f 1 \ / a b 
d^{ c d )-\G 2 )\ c d 

from which it follows, writing G 1 = dG/dr, that 

dX + /dr = (G'/G)X_ + 4GZ = {G'/G, 2G) ■ (X_,2Z). 

But d/dr = fG~ 1 (d/dt) from equation (Ti"4"]) and hence, noting the definition of 
G in (|13p and making use of the Schwarz inequality, one finds 



X+ = (71 - 72, 2/) • 2Z) < [(71 - 72) 2 + 4/ 2 ] 1 / 2 [X 2 + 4Z 2 ] 1 / 2 . 
Combining this result with the above expression for X\ — X^_ then yields 

[x 2 + -4]- 1 / 2 i: + < [(7l _ 72)2 + 4 /2] i/2 j 

which may be integrated to give 

cosh _1 X + /2 < / , 
Jo 



da [(7i ~72) 2 +4/ 2 



Noting L = GX + /2, one finally obtains the Lagrangian inequality 



L < Gcosh 



ds [( 7l - 72 ) 2 + 4/ 



2 , 4f 2il/2 



(29) 
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A sufficient condition for complete positivity in the interaction picture fol- 
lows immediately from (fl2]l. ([27]) , |28|) and ([29j) as 



A < 1, A cosh 



ds[(7i~72) 2 +4/ 2 ] 1 / 2 



< l + A 2 



-2|A-A|. (30) 



It is important to emphasise that this condition can be checked whether or 
not the master equation can be explicitly solved, as it depends only on the 
decoherence matrix 7 in the interaction picture. 

Note that the above condition is 'tight' in the sense that it is in fact necessary 
and sufficient in some cases. For example, when 71 — 72 = / = in the 
interaction picture [6l [8], then G = 1, a = d = 1, and b = c — 0, implying 
that equality holds in (|29| . It is in fact also 'universal', i.e., it is valid for h ^ 
as well as for h = 0, as will be shown below. 



5.4 Generalisations via form-invariance 



While the sufficient condition (j30|) for complete positivity was only derived for 
the interaction picture, with h = 0, it is in fact invariant under rotations about 
the z-axis, and hence may be applied to any master equation with a block- 
diagonal decoherence matrix. 

In particular, for the general case where h is an arbitrary function of time, 
consider a rotation about the z-axis such that h — 0, corresponding to the choice 
a(t) = —(1/2) J Q ds h (see section 4.1). Due to the rotational form-invariance 
of the master equation, it follows that condition (|3H|) must hold with respect to 
the associated decoherence matrix 7, i.e., one has the sufficient condition 



A<1, A cosh 



rf.s[(7i-72) 2 +4/ 2 ] 1/2 



< 1 + A 



,-.2 



-2\A- A| 



for complete positivity. However, using the tranformation equations for 71, 
72, etc. in section 4.1, one finds that all the relevant quantities are rotationally 
invariant. In particular, one has A = A, A = A, W3 — W3, and (71 — 72) 2 +4/ 2 = 
(71 — 72) 2 + 4/ 2 for any choice of a. Hence, condition (|3"D)) is in fact universal. 

The universal form of ([30]) is a fortunate consequence of choosing to work 
in the interaction picture in section 5.3. If, for example, / = instead of 
h = had been assumed, then one would have obtained the simpler evolution 
dX + /dr = (G'/G)X_, leading to the Hamiltonian inequality 



H < G cosh 



ds 1 71 - 72 1 



analogous to the Lagrangian inequality (|29|) . where H replaces L in (|28|) . Thus 
the integrand in equation (|30[) would have been replaced by the quantity I71 —72 1 , 
which is clearly not universal. However, the rotationally invariant form in (|30l) 
can then still be obtained, by considering an arbitrary rotation about the z-axis. 
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Finally, recalling that X\ = X 2 + Z 2 + 4 > 4, it is worth noting that the 
Lagrangian inequality (|29p may be extended to 

1 < L/uj < cosh 

for the average Lagrangian value of the general inverted oscillator equation 
x — uj 2 x = 0. Similarly, the above Hamiltonian inequality may be extended to 

1 < H/u < cosh 

for the average Hamiltonian value of the general oscillator equation x+lo 2 x = 0. 
The upper bounds in these equations can presumably be generalised to include 
an arbitrary function SI, analogous to those for a 2 + b 2 + c 2 + d 2 in [TBJ, leading 
to new sufficient conditions for complete positivity. Minimisation with respect 
to the choice of would then give a 'best' sufficient condition. However, this 
will not be further investigated here. 

5.5 Complete positivity for x?/-damping and zero drift 

The condition that the density matrix remains a positive operator under evolu- 
tion is generally weaker than the requirement of complete positivity, as remarked 
in the introduction. For qubits, positivity is equivalent to the requirement that 
the Bloch vector remains in the unit ball, i.e., |v| < 1. For unital evolution, 
corresponding to the case of a real decoherence matrix 7, and hence w = in 
equation (JBJ, positivity is therefore equivalent to the condition 

M T M < I. 

Consider now a master equation with a decoherence matrix that is both real 
and block-diagonal, i.e., with g = r = s = 0hx equation ([2]). It will further be 
assumed that 73 = 0. These conditions correspond in the Bloch representation 
to zero drift and to the confinement of any damping to the xy-plane. After some 
minor algebra, the condition for positivity reduces in this case to 

A<1, S = x\ +x\ +y\ +yl < 1 + A 2 . (33) 

Noting from ([7|) and (TT2|) that A = A and W3 — 0, this condition is in fact 
equivalent to (|27|) . Hence, complete positivity reduces to the generally weaker 
property of positivity for this particular class of master equations. 

An interesting example of this equivalence is provided by a Bloch-Redfield 
master equation recently studied by Whitney [9], corresponding to a two- level 
system coupled via <7\ to a thermal environment, in the limit of a short memory 
time and high temperature (see equation (22) of [9]). Although this master 
equation cannot be solved explicitly, Whitney has shown via asymptotic analysis 
of |v| that the corresponding evolution satisfies positivity. It follows immediately 
from the above that the evolution must therefore in fact be completely positive, 
despite being unable to construct, for example, an explicit Kraus representation. 



log w + / ds[{Co/uj) 2 + W} 1/2 



(31) 



logw + 


/ ds \(j/u>\ 




Jo 



(32) 
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6 Conclusions 



The main results of this paper are the reduction of a large class of master 
equations to a time-dependent oscillator system and further to a nonlinear 
first-order differential equation, in sections 3 and 4 respectively; the generation 
of new exact solutions using these reductions and rotational form-invariance; 
and the characterisation of complete positivity in terms of an oscillator La- 
grangian/Hamiltonian, leading to a rotationally- invariant sufficient condition 
for complete positivity in section 5. 

Note that further examples of exact solutions, for master equations having 
decoherence matrices as per equation ^ , can of course be generated from any 
exactly solvable oscillator equation. For example, suppose the oscillator equa- 
tion (|15p can be solved for some invertible 'spring constant' function k(r) with 
inverse /c _1 (r), satisfying /c(0) = 1 (which can be ensured by translating and/or 
rescaling r). Making the anzatz / = 0, it follows immediately from (JT3J) that 
f = —h/G = (d/dt)[k^ 1 (G 2 )]. Hence, using the procedure in section 3, one can 
solve the master equation corresponding to the case 

/ = 0, h = -G(d/dt)[k- 1 (G 2 )) 

in equation ([2]), for arbitrary 71 and 72, where G is defined in equation (|13|) . 
Thus any single oscillator solution generates solutions for a large class of mas- 
ter equations. Moreover, one can also obtain solutions for the case / ^ by 
applying the above result to the rotated system defined by / = 0, and using 
form-invariance (this is analogous to the / = example discussed in section 5.4). 
For all such solutions it is then possible to check whether or not the complete 
positivity condition (|27|) is satisfied by the corresponding qubit evolution. 

The further reduction of the master equation problem to a first-order dif- 
ferential equation (12~1]) in section 4.2 provides an interesting subject for further 
study, as it implies that one may also formulate the general time-dependent 
oscillator problem in terms of this equation. For example, it can be shown that 
the particular solution found in section 4.3 corresponds to solving the (inverted) 
oscillator equation (fT5|) for fc(r) = — (1— t/K)~ 2 (for times r < K). More gener- 
ally, a solution of the 'zero-frequency picture' equation (|21|) leads to an implicit 
equation for k{r) which must be solved subject to the condition k(Q) = 1. 

Finally, note for the general case of an arbitary decoherence matrix 7 and 
Hamiltonian H in equation (fT]), the question of complete positivity is invariant 
under unitary transformations, and hence under arbitrary rotations of 7 and 
translations of H (see section 4.1). This implies that the necessary and suffi- 
cient conditions for complete positivity, when formulated in terms of the master 
equation parameters, must be expressible purely in terms of functional of the 
three eigenvalues of 7, with no dependence on H. Now, for the case of a block- 
diagonal decoherence matrix, i.e., of form {2} with r = s = 0, these eigenvalues 
are given by 

A± = ^(71 + 72) ± \ [(71 - lif + 4/ 2 + 4 5 2 ] 1/2 , A 3 = 73. 
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However, the sufficient condition in equation ([50)) can only be written in terms of 
these eigenvalues for the case g = 0. It follows that a stronger sufficient condition 
must exist in general. It is conjectured, comparing the above eigenvalues with 
the forms of equations (TT2")) and (T2SJ), that such a condition may be obtained by 
replacing f 2 by f 2 + g 2 and W3 by w' in equation (|5D|) , where w' is defined by 
the substitution of A+ — A_ for g in equation (|2l)|). 

Acknowledgment: I am grateful to Rob Whitney for motivating my inter- 
est in 'Bloch-diagonal' master equations. 
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